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We suggest a Hamiltonian formulation on a momentum lattice using a physically motivated 
regularization using the Breit-frame which links the maximal parton number to the lattice size. 
This scheme restricts parton momenta to positive values in each spatial direction. This leads to a 
drastic reduction of degrees of freedom compared to a regularization in the rest frame (center at 
zero momentum). We discuss the computation of physical observables like (i) mass spectrum in 
the critical region, (ii) structure and distribution functions, (iii) 5-matrix, (iv) finite temperature 
and finite density thermodynamics in the Breit-frame regularization. For the scalar <f)3^i theory we 
present numerical results for the mass spectrum in the critical region. We observe scaling behavior 
for the mass of the ground state and for some higher lying states. We compare our results with 
renormalization group results by Liischer and Weisz. Using the Breit-frame, we calculate for QCD 
the relation between the FF'"' tensor, structure functions (polarized and unpolarized) and quark 
distribution functions. We use the improved parton-model with a scale dependence and take into 
account a non-zero parton mass. In the Bjorken limes we find the standard relations between Fi, 
F2, gi and the quark distribution functions. We discuss the role of helicity. We present numerical 
results for parton distribution functions in the scalar model. For the c^^-model we find no bound 
state with internal parton structure. For the c^^-model we find a distribution function with parton 
structure similar to Altarelli-Parisi behavior of QCD. 

PACS-index: 13.85.-t, ll.lO.Ef 



I. INTRODUCTION 

The standard model of strong interaction physics (QCD) has been confirmed very successfully by comparison 
between experiment with perturbative and non-perturbative (mostly lattice) calculations. There are a number of 
observables, which need to be computed non-perturbatively. In some cases non-perturbative computational progress 
seems very hard to come by. Let us mention the following examples: (a) S'-matrix for hadron-hadron scattering, (b) 
Structure functions of the proton, in particular at small values of the Bjorken variable xb. (c) Excited states in the 
hadron mass spectrum, (d) Finite density thermodynamics of hadronic matter. 

In deep-inelastic lepton-hadron scattering (DIS) one is interested in structure functions and their interpretation in 
terms of distribution functions. In order to define distribution functions, one has to specify a reference frame. Possible 
choices are the rest frame, the infinite momentum frame, light-cone coordinates, or the Breit frame. Conventionally, 
the infinite momentum frame and light-cone coordinates have been mostly used. The Breit-frame (characterized by 
go = 0) is a distinguished frame in the sense that Q = \/—q^, the momentum of the exchanged photon is the resolution 
ability of the photon to resolve the proton structure. This does not hold in any other frame with go 7^ 0. 

The Breit-frame is not only conceptionally attractive, but we suggest here that is is useful also for non-perturbative 
numerical calculations. In Ref. [1] we have introduced a new regularization scheme for a lattice Hamiltonian, based on 
the Breit-frame. Its construction is guided by the kinematical variables which play a role in deep-inelastic scattering 
(DIS). We have introduced a momentum lattice based on a the Breit-frame. It is centered around the proton 
momentum (in the case of DIS proton scattering). The scheme restricts parton momenta to positive values in each 
spatial direction and links the maximal parton number to the lattice size. This leads to a drastic reduction of degrees of 
freedom compared to a regularization in the rest frame with the center at zero momentum. In Ref. [1] we have applied 
the scheme to the scalar (/)|_|_]^ -model. We have computed the mass spectrum and extracted physics close to the critical 
line (second order phase transition). We found very good agreement with the predictions of the renormalization group 
by Liischer and Weisz [2]. To our knowledge, in Ref. [1] critical behavior of a 3 -|- 1 dimensional field theory has been 
extracted for the first time from a Hamiltonian formulation. 

The successful working of the method with respect to the relative small numerical effort (diagonalization of matrices 
in the order of 50) makes us cautiously optimistic that other physical observables or other models could be treated as 
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well. In this work we want to elaborate on these ideas. In sect. 2 we present the Hamiltonian formulation in the Breit- 
frame regularization. We explain physical and mathematical reasons for the working of the method. The calculation 
of the mass spectrum for the 4''^^i model with numerical results near the critical line are presented in sect. 3. In sect. 4 
we discuss the structure functions and distribution functions in the Breit-frame. For QCD we compute analytically 
the relation between the hadronic tensor W^'^ , the structure functions Fi, F2, gi, §2 and the quark distribution 
functions. We present numerical results for the parton distribution function for the scalar models <f)'^^i and (fis^i. The 
usefulness of this method eventually depends on its potential in numerical calculations of gauge theories. Thus the 
Breit-frame regularization for lattice gauge theories is given in sect. 5. In sect. 6, we discuss advantages of the Breit- 
frame regularization for the purpose to compute the S'-matrix of a scattering reaction from a Hamiltonian lattice 
formulation. Finally, the computation of thermodynamical observables at finite temperature and finite density from 
the Breit-frame is the topic of sect. 7. A summary is given in sect. 8. 



Considering non-perturbative methods in many-body physics, statistical mechanics and field theory, most successful 
techniques are sum-rule techniques and lattice field theory in the Lagrangian formulation using Monte-Carlo methods 
to compute functional integrals. Although there is a Hamiltonian formulation of lattice field theory, i.e., the Kogut- 
Susskind Hamiltonian in lattice gauge theory [3], Hamiltonian methods have not been mainstream in the domain of 
non-perturbative methods. One of the basic reasons was that in order to treat adequately the physical degrees of 
freedom, very many virtual particles have to be taken into account. As a function of particle number the dimension of 
Hilbert space increases exponentially. Nevertheless (maybe due to shortcomings or slow progress in Lagrangian lattice 
field theory), over recent years several workers have explored Hamiltonian methods, by trying to work with effective 
Hamiltonians with a small number of degrees of freedom. Examples are the work by Liischer [4] and van Baal [5], the 
exp[5'] method coming from nuclear physics [6], applications of the Kogut-Susskind Hamiltonian to compute glueball 
masses and string tension in QCD [7], the Hamiltonian approach in light-cone quantization [8], and quite recently a 
Hamiltonian renormalization group approach [9]. These approaches have employed quite different strategies to cope 
with the problem of a large number of degrees of freedom: E.g., Wilson and co-workers have pursued the idea of the 
renormalization group, i.e., thinning out degrees of freedom and constructing a new (renormalized) Hamiltonian with 
a sufficiently small number of effective degrees of freedom. Liischer [4] and van Baal [5] have discovered that much 
physics of the low-lying QC-D-spectrum, at least for small lattices can be described by zero-momentum dynamics 
plus a suitable treatment of the remaining degrees of freedom. The idea of the exp[5'] method [6] is that the linked 
cluster structure underlying a ground-state in a many-body theory can be generated by a suitable operator S, and 
automatically guarantees the correct infinite volume singularity of the ground state energy. In the applications of 
the Kogut-Susskind Hamiltonian to QCD [3], several groups have developed clever ways to take into account the 
high number of plaquettes and closed loop variables, e.g., via the t-expansion method by Horn and co-workers [10]. 
Finally, an advantage of the regularized (discretized) light-cone Hamiltonian method is that light-cone momenta p'^ 
of all partons are positive and add up. Thus a total light-cone momentum P'^ drastically limits the number of degrees 
of freedom. However, this does not hold for the perpendicular momentum 

Here we pursue the following alternative. Let us consider as example the scalar model given by the Hamiltonian 



where mo and go are the bare mass and coupling constant, respectively. We construct the corresponding Hamiltonian 
in momentum space 



II. FORMALISM 



A. Hamiltonian formulation 
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We did not normal order, but we have subtracted the vacuum energy. We have written the Hamiltonian in discretized 
form by introducing a lattice in momentum space with a momentum resolution Ak and a momentum cut-off A. 
Conventionally one would choose a regular lattice, being symmetrical with respect to momentom zero (rest-frame). 
This would be suitable to compute the vacuum state with the quantum number P = 0. Our alternative is: We choose 
the same regularization, however, we retain only those lattice momenta which correspond to fast moving partons 
going in the same direction as the hadron (proton). This is called Breit-frame regularization in the following. As will 
be shown below this yields a drastic reduction of the effective degrees of freedom compared to the rest frame. 
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B. Breit-frame regularization 

The most important experiment in order to probe the structure of hadrons is deep inelastic scattering (DIS). 

Its simplest form involves inclusive scattering of an unpolarized lepton off a hadronic target. Let us recall some 

basic notations [11]. The hadron in its ground state interacts with the probing lepton by the exchange of a virtual 

photon (or neutrino). The hadron (proton) carries momentum P before the collision and goes over to a hadronic 

state X with momentum Px ■ The electron has correspondingly momenta k and k' . The exchanged photon carries 

momentum q = k — k' . One defines = —q^. In Feynman's parton model it is assumed that the proton is made up 

of constituents, i.e., the partons. They are weakly bound, i.e., the binding energy is small compared to the resolution 

ability Q := ^/—qjlq^ of the probing photon. Let p denote a parton momentum. Conventionally, one introduces the 

q2 ^ p p{L) 

Bjorken variable xb '■= = = — ttt. The superscript L denotes the longitudial direction, i.e. the direction 

^PfiH^ q ■ P P^ ' 

of P. The second equation results from the impulse approximation, i.e. the partons are on the mass shell. The last 
identity holds in the Breit-frame. The Brett frame of the hadron is defined by the requirements that the photon 
energy go be zero and that the photon momentum q be antiparallel to the hadron momentum P . In this frame, the 
longitudinal component of the parton momentum obeys = '5/2. The rationale for this particular choice of frame is 
that QCD structure functions F[xb,Q) can be interpreted as a linear combination of parton momentum distribution 
functions f{xB, Q), which have a more intuitive interpretation. This relation holds for leading twist (higher twists are 
suppressed for large Q^). Structure functions are another way of expressing scattering cross sections. The distribution 
function of a parton counts the number of those partons with a given momentum fraction a;^ in the proton. For a 
precise definition see Ref. [11]. 

If the hadron is in its ground state, then the longitudinal momentum of the parton can neither be negative nor 
can it be greater than P^. 

0<pW<pW, (3) 

Thus follows the well-known constraint < xb < I- Eq.(3) can be viewed as a regularization of the longitudinal 
parton momenta. However, it does not restrict the transverse components of parton momenta. In the Bjorken limit 
[Q^ oo, P ■ q ^ oo, Xb = const.) combined with the Breit-frame the hadron is a fast moving object. In momentum 
space, an object which is spherical in the rest frame becomes prolate in a fast moving frame. Hence it is physically 
justified to restrict the transverse parton momenta to a finite region of ellipsoid (prolate) shape. In particular, we 
have chosen a sphere centered at the mid-point of the interval [0, P'^^^ (any prolate ellipsoid lies within this sphere), 
given by 

ip-P/2f <iP/2f. (4) 

One should note that this constraint also follows directly from < xb < I and going into the "parton Breit-frame" 
(defined by g*^ = and being anti-parallel to the parton momentum p) where xb takes the form xb = -p^- 

Because we are working in the Hamiltonian approach we need to define a basis in Hilbert space. We construct the 
Hilbert space as a Fock space of free particles and select (parton) momenta p from a bounded domain corresponding 
to DIS as given by Eq.(4). This is an assumption based on the physical intuition that the experimentally observable 



3 



parton momenta are those which dominate the quantum dynamics. This assumption has been tested by computing 
critical behaviour of renormalized masses and a good agreement with analytical scaling behaviour has been observed 
(see below). 

Now we introduce a momentum lattice regularization: In order to have a practically convenient lattice we further 
constrain the parton momenta from Eq.(4), namely by selecting a regular cube centered at _P/2 and located inside 
the ball given by Eq.(4). I.e., the parton momenta]? are restricted to the domain 

\P\ 

< Pi < A = --j= for i = X , y, z . (5) 

Inside this domain we define lattice momenta p := nAp where n is an integer vector and Ap is the momentum 
lattice resolution. One notices that all lattice momenta are non negative. Contrary to regularization in the rest 
frame which does not limit the particle number, our approach has the following important property: The effective 
Hilbert space is built from the Fock states | (at )"i • • • (at | >, with the conditions that the total momentum 

be Tiiki + • • • + niqkiq = P , which is located on the surface of the Breit domain, and that each parton momentum ki 
be inside the Breit domain, given by Eq.(5). Thus the regularized Hamiltonian is given by Eq.(2), restricted to the 
effective Hilbert space. 



C. Reasons for reduction of number of effective degrees of freedom 

Why does this regularization scheme lead to a Hamiltonian with a small number of effective degrees of freedom? 
Firstly, for any given state from the effective Hilbert space, the Fock space particle numbers are bounded, if one 
considers only non-zero parton momenta. This follows from niki + • • • + n^qkiq = P. Thus an upper bound on the 
\p\ 

total particle number is This does, however, not limit the zero-mode particle number. The zero-mode has to be 
taken into account explicitely. The zero-mode only determines the vacuum expectation value of the field < 4> >. In 
this work we only consider the symmetric phase of the model < (f) >= 0. In the model the vacuum expectation 
value < (/) > is an order parameter for symmetry breaking and thus the field has fiuctuation zero in the infinite volume 
limit. In this limit it becomes a classical variable. Thus for sufficiently large volume it is justified to set the zero-mode 
to zero. Besides, in models where the zero-mode can not be dropped, the zero-mode describes only one degree of 
freedom, which can be treated like a quantum mechanical oscillator. In summary, the ultraviolet cutoff A given by 
Eq.(5) implies a total particle number cutoff and thus drastically reduces the dimension of the Hilbert space. 

Secondly, if one wishes to compute the mass spectrum of a physical particle, but does not want to compute the 
vacuum, one has the freedom to choose a reference frame boosted to a momentum P z/z 0. As is well known from 
many-body theory and the exp[5'] method, the vacuum state energy has a volume divergence, but the energy of a 
physical particle state does not have such a divergence. Thus choosing a sector with P z/z excludes the vacuum 
state, but may eventually allow to compute more easily the mass of a physical particle, compared to a computation 
in the rest frame where the vacuum is present. 



III. MASS SPECTRUM AND CRITICAL BEHAVIOUR OF <^|+i THEORY 

Firstly, we need to convince ourselves that the method allows the correct computation of physical observables. 
We have chosen the scalar 4''^^i theory because it is a quite well understood theory which has a second order phase 
transition, allowing to test our method near a critical point. The Hamiltonian of the (j>^ theory is given by Eq.(2), 
constrained by the Breit condition (5). It is expressed in terms of free field creation and annihilation operators 
corresponding to the lattice momenta in the Breit-frame. Because the Hamiltonian and the momentum operators 
commute, we compute the energy spectrum E„ in a Hilbert space sector of given momentum P. Since we are not in 

the rest frame we use the mass-shell condition M„ := — P^ in order to obtain the physical mass spectrum. 

It is known [2] that the critical line between the symmetric and the broken phase lies entirely in the region where 
the bare parton mass squared mg is negative. Hence we cannot build the Fock-space in terms of partons with those 
masses. As a remedy we have split the bare mass squared mg = ni^j^ -|- into a positive kinetic part ni^j^ and an 
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interaction part The Fock states are built from positive bare masses rrikin- In numerical calculations close to 

the critical point shown in Fig.[l] we have chosen, for simplicity, a small positive value. We found that the lower lying 
physical mass spectrum is not very sensitive to the value of ffikin 

(this is not the case for higher lying masses). A 
better choice of nikin in our view would be be to take the renormalized mass niR. Although niR is unknown initially, 
it can be computed by making an initial guess and then iteratively improving the answer. 

We have diagonalized the Hamiltonian on two lattices: A/Ap = 3 and A/Ap = 4. This would correspond to 
symmetric lattices [—A, +A] of size 7^ and 9^ nodes, respectively. This results in a very small Hilbert space of only 
6 and 21 states, respectively. In order to compare our results to those of Liischer and Weisz [2] we express the bare 
parameters mo and go in terms of the parameters A and k: 

m2 = (1 - 2A)/«: - 8, g, = &\. (6) 

Fig.[l] displays the renormalized mass niR versus k. One observes that our results computed on very small lattices 
are quite close to the results of Liischer and Weisz [2]. Masses M computed on the lattice must obey a < l/M < L, 
where L is the length of the lattice and a denotes the lattice spacing of a space-time lattice. It is related to A by 
A = It can be shown from perturbation theory [12,2] that the physical masses close to the critical point obey the 
following scaling law 

M ~ Crl/2|/r^r^l/^ (7) 

where r := 1 — k / Kcrit and C is a constant (integration constant of renormalization group equations) . Since the results 
of Ref. [2] are based on the solution of the renormalization group equations, this scaling law fits their results. One 
should note, however, that two different regularizations (this work and that of Ref. [2]) in general correspond to two 
different critical lines corresponding in general to different values of Kcrit- In Tab. [1] we have displayed our results 
for the critical points Kcrit as a function of A and compared our results with those of Ref. [2]. Again, our results are 
very close to those of Liischer and Weisz. These results cover a domain of the bare parameter space extending quite 
far away from the Gaussian fixed point at k = 1/8, A = 0. 

Another way to test continuum physics is to look at the mass ratios M„/Mi from the spectrum on the lattice 
and check if they become independent of the cutoff A or else independent of the coupling constant (j'o(A) (i.e. they 
scale). These mass ratios M„/Mi are shown in Fig. [2]. As can be seen, for a number of states M„/Mi const in 
a wide range of k- values, i.e. they scale. However, for some states M„/Mi diverges, i.e., there is no scaling. The 
physical reason behind this is the following: The 4't+i niodel describes a gas of partons repelling each other [2]. The 
spectrum of Fig. [2] shows states dominated by the l-,2-,3-,4- particle Fock space sectors plus a spectrum of excited 
(scattering) states. The picture of repulsive two-particle-exchange force is confirmed by observation that the mass 
of the lowest-lying n-body state is larger than n-times the mass of the one-body state. The states which scale are 
just those lowest-lying n-body states. The higher-lying part of the spectrum consists of states with more nodes in 
the wave-function than lattice points, having also a wider range and contributions from higher Fock-state sectors. 
Because in the calculation corresponding to Fig. [2], the parameters Ap, A and the parton number cutoff are all kept 
fixed, we cannot properly describe these higher-lying states. Consequently, they do not show scaling. When we go to 
bigger lattices (Ap 0) we then observe (not displayed here) more states which show scaling. 

IV. STRUCTURE FUNCTIONS 

A. Why structure functions in the Breit-frame? 

Hadron structure is probed by deep inelastic scattering (DIS). Over recent years a great deal of experimental 
data has been gathered from high energy collider experiments. While perturbative quantum chromodynamics (QCD) 
describes successfully the large dependence of DIS structure functions, it cannot predict the correct dependence 
on the Bjorken variable xb- Thus much effort has been devoted to compute quark or gluon distribution functions 
and proton structure functions from QCD with non-perturhattve methods. E.g., Martinelli et al. [13] have computed 
the first two moments of the pion structure function via Monte Carlo lattice simulations. Schierholz and co-workers 
[14] have recently computed moments of proton and neutron structure functions. These calculations are notoriously 
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difficult. A particular problem is the determination of small a;B-behavior from a few moments (for the present status 
of lattice calculations of structure functions see Ref. [15]). This situation calls for alternative techniques. 

Let us briefly outline the reasons for the choice of our method: (i) Structure functions are computed from wave 
functions. Wave functions are defined in Minkowski space. The Hamiltonian approach offers the advantage of allowing 
direct computation of Minkowsky space observables. E.g., scattering wave functions for glueball-like states in compact 
QED2+1 have been computed in a Hamiltonian formulation on a momentum lattice [16] (for a review of Hamiltonian 
lattice methods see [17,7,8]). (ii) The usefulness of a momentum lattice to compute physics close to a critical point 
has been demonstrated in Ref. [19-24]. (iii) The reason for our choice of the Breit frame has been explained above. 
However, Hamiltonian methods are known to lead to numerical problems because of the huge number of degrees of 
freedom involved [25]. To the best of the authors' knowledge nobody has succeeded before to observe scaling behaviour 
indicating continuum physics in a (3+l)-dimensional Hamiltonian lattice formulation. 

The Breit frame has a distinct property: Only in this frame the photon momentum transfer Q can be interpreted 
as resolution ability of the photon. The quark and gluon distribution functions of a proton or a neutron which are 
measured by DIS show a peak for small xb even for a moderate resolution Q [26]. This indicates a huge number of 
partons in the proton, because a system of n identical partons would be peaked at xb = l/« for symmetry reasons. 
The physical reason is that the strong forces which bind the proton can easily create gluons or quark-antiquark 
pairs. Contrary to a typical non-relativistic problem, particle number is not conserved. Consequently, because of 
the enormous number of degrees of freedom which are usually associated with a relativistic many-body system, it is 
almost impossible to calculate quark or gluon distribution functions or mass spectra in a Hamiltonian QCD approach 
[27]. Our regularization, given by Eq.(5), however, enables us to treat a large number of partons with a reasonable 
numerical effort. 

The Q dependence of the distribution functions is also a many-particle effect. An intuitive explanation for this 
dependence is that more partons can be seen inside the proton, if the resolution Q is increased. Partons, however, 
which are heavy with respect to the forces between them, can be described in a simple constituent model because 
many-particle effects are negligible. Hence, their distribution functions are neither peaked at xb = nor do they 
depend considerably on the resolution Q. Examples are heavy quarkonia, electromagnetically bound particles (such 
as atoms) or the 4''^^i theory which we are investigating below. The renormalized coupling constant of 4''^^i theory 
is weak everywhere in the critical region and the forces between "partons" are even repulsive [2]. 



B. Relation between structure functions and distribution functions in the Breit-frame for QCD: unpolarized 

structure functions 

In this section we compute analytically the relation between hadronic tensor, structure functions and distribution 
functions. Because we work in a fast moving frame and not in the infinite momentum frame, we can take into account 
explicitly a non-zero parton mass. The cross section for deep inelastic lepton-hadron (electron-proton) scattering has 
the form (following Jaffe's notation [28]) 



d'aocl'^'^W.,, 



(8) 



where /^j^ denotes the leptonic tensor and Wf^, stands for the hadronic tensor. The hadronic tensor can be split into 
a symmetric part, which corresponds to unpolarized structure functions and an anti-symmetric part, corresponding 
to the polarized structure functions. The symmetric part can be parametrized in terms of the structure functions Fi 
and F2, 



*^ *^ sym 



Fi 



El 

1 

V 



(9) 



where q^ is the photon momentum, is the proton momentum and v = q ■ P . Now we choose as reference frame 
the Breit-frame: In the Breit-frame the proton momentum is P^ = (£',0,0,-P3) with = P^ + M^, M being the 
proton rest mass. The photon momentum is qf = (0, 0, 0, —Q), Q is defined to be q^ = —Q^. As a result we find that 
all components of W^^^ vanish, except for 
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J7'2 

W°° = -F-i H F-> 

* * sym P Q 

VI/ii = W^^ = F^ flO) 

* * sym * * sym -L V / 

The hadronic tensor is defined [28] as 

iTTW^"" = ^(27r)4(5(P + q - Px) < PS \ J^(0) \ X X X \ .r(0) \ PS >, (11) 

X 

where X denotes the unobserved fragments of the proton, P is the proton momentum. We have normalized the proton 
state to < P' \ P >= 2E{2it)^5^{P' - P). S \s the proton spin (Pauli-Lubanski vector) normalized to S'^ = -M'^ . 
Jt'(x) = ip{x)'jf'ip{x) denotes the fermionic (quark) current. The hadronic tensor can be expressed as a current 
commutator, 

iTTW''^ = j d^yeM-n ■y]<P\ [J'iy), rm \p>- (12) 

In deep inelastic scattering it is customary to use the impulse approximation. The partons lie on the mass shell. Thus 
one can expand the field il>{x) 

= E/ (^(2^W)"'^' [u.{k)e-^'-h,[k) + v,[k)e^^-d\[k)\ . (13) 

The spinors are normalized to UgUgt = 2m5g gt, VgVgt = —2m5g gt, where m is the parton rest mass. The parton 
spin is normalized to = — m^. The creation and annihilation operators obey [&s(fc), &|,(fc')]_|_ = 6g gi6[k — k') and 
[dg{k),dl,{k')]+=5g^gi5{k-k'). 

In the computation of the matrix element of the current commutator the following leptonic tensor occurs: 
lannnik, s, k' , s') = u{k , s)-f''u{k', s')u{k', s')-f''u{k, s) 

= k^k'" + k'^k" + g^^im' -k-k')- ime"""!^ {k - k')„{s + s')/3, (14) 

and summing over s' yields 

Cau(^, k', s) = 2 [k^k'" + fc'^r + g^^im' -k-k')- ime'"'"'^{k - fc')aS/3] , (15) 

being the standard result [28]. Due to the current commutator, there are four fermion fields involved, which gives 
16 combinations of fermion and anti-fermion creation and annihilation oprators. A straightforward but lengthy 
calculation gives the following term (other terms see below) 

w,Z,,, = j f ^3fe '^(g° + fc° - n j-i^^^^j,^^ ^ ,)n5(p^ fc). (16) 

J uj[k)uj[k + q) ^ 

Here we have switched to the following normalization of the proton state < P' \ P >= 5^[P' — P). In the calculation 
occurs the the matrix element < PS \ h\(k)hsi(k') \ PS > which allows to split off S[k — k') due to total three- 
momentum conservation as well as to split off Sg^s' due to conservation of spin quantum numbers in the helicity basis. 
Thus we have defined Il''g[PS, k) by 

< PS I bl{k')b,{k)PS >= U'l{PS,k) S{k-k')S,y. (17) 

which is the the expectation value in the proton state of the quark number operator corresponding to momentum k 
and spin s. Note that 11* has the same dimension as < PS \ PS >. Now we go into the Breit-frame. In particular, 
we employ the Breit-condition, Eq.(4) and (5). In the Breit-frame one has q° = 0, moreover k lies on the mass shell, 
(fc'^)2 = -\- , due to the impulse approximation. Also k' = k -\- q \s on-shell. Thus we obtain 

Kn.n = f / d^^^ E %^^i|^nt(P, k,), (18) 
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where we have defined kq = {k±,Q/2,) and k'q = {k±, —Q/2), (parallel and perpendicular denotes the orientation 
of components with respect to the space-component of the proton momentum). 



A second term, which contributes to the W^"^ tensor is the following 



* J , Quj(k'o) 



Ilf{PS,k) is defined in analogy to Eq.(17), but for the anti-quark number operator. The leptonic tensor l^yyy 
corresponding to the t;-spinor is defined by 

l^:,y{k,s,k',s') = v{k,sh''v{k',s')v{k',s')Yv{k,s), (20) 

yielding, after summation over the spin s' 

l^yyy{k, k', s) = 2 [k^k'" + k' " k" + g {m^ -k-k') + ime'"'"'^{k - fc')aS/3] • (21) 

All other terms give vanishing contributions to W^"^ due to the fact that all parton momenta lie in the Breit-sphere. 

Thus the symmetric part of W^'^ , corresponding to the unpolarized structure functions yields the following result 
in the Breit-frame. The only non-zero components are those with jj, = v = Q,l,2. 

yV^'m = f / E k(P, kQ) + nf (P, kQ)] , (22) 

where l^npoi denotes the symmetric part of l^yyy and l^yyy 

Cnpoiik, k') = 2 [k^k"^ + k'^k'^ + g^^m' - k ■ k')] . (23) 

In particular one has 

Cnpol{kQ,h'Q)=A[m' + {k^fl 

ll\poi{kQ,k'Q)=A{k,f + Q\ 

lllpo!ikQ,k'Q)=4{k,f + Q''. (24) 

Let us now consider the Bjorken limes of those expressions. The Bjorken limes is defined by Q — cxd and xb = const. 
In the Breit-frame, this implies for the proton momentum Ps that Ps ^ co and Q / P3 = 2xb = const. Thus we 
compute the kinematical factors E = \/ + P^ ~ Q/2xb and 



This implies 



Q-^^O. (25) 



W^ym ^ 0. (26) 
On the other hand, Eq.(lO) yields in the Bjorken limes 

W°° = -Fi + ■^F2 ^ -Fi + -^F2. (27) 
P3Q 2xB 

Thus the last two equations imply in the Bjorken limes that the Callan-Gross relation holds 

2xbFi = F2, (28) 
which is the standard result as in the parton model. Similarly we compute 
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Thus we obtain the following result for W'^y^, ^sym ^^'^ ^svm' Bjorken limes 



(29) 



14^°° =0 

sym ' 

p 

Pyll = Py22 _3 

sym sym 



j d'k^ J2 + n^(p, kQ)] . (30) 



Our regularization scheme allows to compute directly the parton distribution function. We define 

f(p3, Ps, l^)= j Cpk^ J2 kg) + n^(P, kg)] , (31) 

where ps is the longitudinal parton momentum, Ps is the longitudinal proton momentum and is a fixed but 
arbitrary scale parameter with dimension of mass (e.g., Aqcd)- f{P3, P3, Ij) is the probabilty of finding a parton with 
longitudinal momentum in a boundstate (proton) with longitudinal momentum ^3, where momenta are measured 
in terms of the scale ji. The W^'^ tensor is dimensionless. Thus we have 

W^^rr^ = y/- (32) 

Thus G{p3, P3, jj) = P3f{p3, P3, jj) is a dimensionless function which scales 

G(Aj93,AP3,A//) = G(j93,i^3,/«). (33) 

We have shown in the Bjorken limes that W'fym = Pi) thus 

Fi{x,Q,l^) = YfiPs^Ps^l^)- (34) 
Expressing p3 and P3 in terms of Q and x yields P3 = Qj'l and P3 = Qj'lx and hence 

^^i(a;,0,//) = £/(0/2,0/2a;,//). (35) 

Making a scale transformation by multiplying all variables of dimension mass by A, where A is chosen to obey Ai^s = 1, 
yields 

Fi{x,Q,^) = ^f{x, l,2x^/Q), 

F'2{x,Q, n) = xf{x,l,2x^lQ). (36) 

Here f[x, l,2xiJ,/Q) denotes the probabilty of finding a parton with longitudinal momentum fraction x from total 
longitudinal momentum = 1, where the scale is given by 2xiJ,/Q. Note that in Eq.(30) the sum runs over all spin 
values. For a spin 1/2 parton this is equivalent to a sum over the helicity quantum numbers + and — . If we take 
into account e,-, the electric charge of quark of fiavor i relative to the charge of the electron, and redefine Eq.(31) by 
/ = / £k_L J2s and / = / £k_L J2s nf we obtain from Eq.(36) 

Fi{x, Q) = ^e^l [f!l\x, 1, 2xfi/Q) + f^_!\x, 1, 2xfi/Q) + f^^^x, 1, 2xfi/Q) + f^_!\x, 1, 2xfi/Q)] 

i 

F2{x, Q) = ^e^a; [fl\x, 1, 2xii/Q) + f1\x, 1, 2xii/Q) + fl\x, 1, 2xii/Q) + fi\x, 1, 2xii/Q)\ . (37) 

i 

The standard expression from the parton model [11,28] is given by 
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i 

F,(x) = Y.e}x [ql\x) + q^^\x) + c^{> (x) + -q^^\x)\ . (38) 

i 

In the Bjorken limes our result, Eqs.(37), agrees with the standard result, Eq.(38). As can be seen, the quark distri- 
bution function q[x) occuring in Eq.(38) does not have any Q dependence. It corresponds to the naive parton model, 
which has no Q dependence. However, perturbative QCD introduces a Q dependence via logarithmic corrections 
(violation of scaling). Thus one arrives at q{x,Q,iJ,) which is a quark distribution function from a "renormalization 
group improved parton model". q{x,Q,iJ,) is interpreted as the probability to find a parton with momentum fraction 
a; in a hadron with momentum = cxd, where the resolution (by the photon) Q is finite, and momenta are measured 
in terms of a mass scale Note that our distribution function f{x, l,2xiJ,/Q) has a different interpretation than 
q{x,Q, jj): f{x, l,2xiJ,/Q) corresponds to the Breit frame where the hadron moves fast but with finite momentum, 
while q{x,Q,iJ,) corresponds to the infinite momentum frame. However, in the Bjorken limes both coincide. 



C. Relation between structure functions and distribution functions in the Breit-frame for QCD: polarized 

structure functions 

The anti-symmetric part of the W^'^ tensor describing the spin dependent part can be parametrized in terms of the 
spin structure functions gi and g2 [28] 



-ie^""'Pq„ 



Sp , . q ■ SPp 
— (91 +92) 5-^ 



92 



(39) 



Here S denotes the proton spin, with P ■ S = and = —M^. The proton spin can be polarized in two ways: 5 || -P 
(longitudinal, helicity) or 5 -L -P (transversal). In order to extract both spin structure functions from the tensor W^^^ 
one needs both polarizations. Let us consider first longitudinal polarization. Then we find in the Breit-frame 



W, 



12 



91 



92 



W. 



21 



-w, 



12 



(40) 



All other elements of W]^]! vanish. Now let us consider transversal polarization. Then we find in the Breit-frame 

M 



'^:t = -i^{9i+9'2) 



20 



-w, 



02 



(41) 



All other elements of W]^]! vanish again. 



In the following let us consider the case where the proton as well as the partons are polarized longitudinally (helicity). 
We want to compute in the Breit-frame the 12-component of the hadronic tensor W^^ . We obtain 



^pol i^Q' ^'q' 



Q. 



(42) 



where we have defined 



'po; (^Q ) ) s) — ^auUui^Q J ^qj ^) ■ 



(43) 



In the Breit-frame the anti-symmetric part of the leptonic tensor l^uuui^j *)) Eq.(15), and oil'^yyy[k, k' , s), Eq.(20), 
are given by 



^auau i^j k , s) — i2Qso, 



(44) 
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For the parton spin s holds s ■ k = 0, s"^ = —m^. In the helicity basis, one has s \\ k. thus one defines the parton 
helicity h by 



In the Bjorken limes one obtains 



and hence 



* as 2 



h. = —jr- (45) 
ui \ k \ 

;i2 

— S^^i2h,. (46) 



Eq.(40) implies in the Bjorken limes 



W'J^igi. (48) 
From this and Eq.(47), after doing the same scale change as in the unpolarized case, we arrive at 

Q) = IJ2 [f'+f+i^' 1' 2a;///0) + h.f^_!\x, 1, 2xfi/Q) 

i 

+ h+ff{x,l,2xii/Q) + h_fl\x,l,2xii/Q)\. (49) 
This is in agreement with the standard result of the parton model [28] 

5l(^) = ^E^' [h+ql\x) + h.q^l\x) + h+-ql\x) + h.-q^l\x)\ . (50) 

i 

In summary of this section, we have computed analytically in Breit-frame regularization the relation between 
hadronic tensor, structure functions and parton distribution functions. The main results are given in Eqs. (22,23) 
and Eqs. (42, 43). The results are "renormalization group improved" compared to the naive parton model, taking into 
account parton mass m and scale parameter ji. 



D. Numerical results for distribution functions of the scalar model 



In this section we want to show how distribution functions can be computed numerically with the Hamiltonian 
approach in the Breit-frame regularization. We apply the method to the scalar model in 3 + 1 dimensions. This model 
has been extensively studied and represents, for finite cut-off, a non-trivial effective theory. We compute distribution 
functions for the and (f)^ model. Let us firstly consider the (t)'k,j^i model. The first excited state consists of only 
one parton (three-particle contributions have been found to be extremely small, i.e. within the error margin). This 
is because our particular choice of regularization (Breit-frame), which makes all parton momenta positive, implies 
for the Hamiltonian that those terms are dominant which conserve particle number. Consequently, the distribution 
function is peaked at xb = 1, as shown in Fig. [3]. We did not find a noticeable dependence on the resolution Q, i.e., 
many-particle effects are absent. We have not put more effort to obtain a finer xb resolution, because this state does 
not display the interesting structure of a bound state. We have observed also that higher excited states display a 
dominant 2-, 3-, 4- particle content (with very small mixtures between different sectors). One should note, however, 
that the simplicity of the first excited state is due to the fact that the positivity of the longitudinal parton momenta 
prevents the creation of partons directly from the vacuum. Had we worked in the rest frame (_P = 0), the "valence 
parton" of the first excited state would be surrounded by a large cloud of partons with opposite momenta p and even 
the vacuum, lying in the _P = sector, would be made up of such a cloud. 

Now let us consider the (f)^ model. We are aware that this theory is mathematically not well defined: It suffers from 
an unstable vacuum since it is unbounded from below. On the other hand, it is well defined in perturbation theory. 
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It is interesting to note that the (f)^ model in _D = 6 dimensions has the property of asymptotic freedom [18], opening 
the possibility to study perturbatively scaling violations. In _D = 4 dimensions asymptotic freedom is known to exist 
only for non-Abelian gauge theories. This (t)%j^i model serves here to illustrates how Altarelli-Parisi like behaviour 
and a sharp forward peak for small xb can be obtained in our method. In hadron physics one is ultimately interested 
in distribution functions of the proton which is a bound state. Its distribution functions reflect the strong variation 
of parton number of the proton and they look quite different from what we have seen above for the (j)"^ model. Thus 
in order to illustrate the capability of our method to treat such strong particle creation effects with a reasonable 
numerical effort, we have chosen the scalar (f)^ model, which yields a bound state. We are aware that this theory is 
mathematically not well deflned: It suffers from an unstable vacuum since it is unbounded from below. This unstable 
vacuum prevents a meaningful calculation of ground state masses which are needed to specify renormalization group 
trajectories and hence the exact relation between the resolution Q and the bare coupling constant go. Therefore we 
can only compute the distribution function f{xB,go) where go is the (j>^ coupling constant. We have computed the 
distribution function in l-,2- and 3 space dimensions. For a given parton number cutoff these curves look very much 
alike. We present the result corresponding to a calculation in one space dimension (Fig. [4]) with A/Ap = 11 which 
implies that up to 11 partons can be present in the Fock space. When the coupling go increases we see that the 
distribution function develops a peak at momentum fraction xb = 0- This is because increasing the coupling means 
that more partons are produced which share the total momentum fraction. Decreasing the parton masses produces 
the same effect. Without a suitable regularization, the numerical effort to describe a system with a large particle 
number would be drastically higher, even for the (j>^ theory. 

V. APPLICATION TO GAUGE THEORIES 

Because the physically most important models are gauge theories we want to discuss the treatment of gauge theories 
in the Hamiltonian formulation with Breit-frame regularization. In the previous sections we have given arguments 
and numerical results showing the usefulness of a momentum lattice regularization in connection with the Breit- 
frame. The usefulness of a momentum lattice corresponding to the rest frame has been previously investigated and 
demonstrated by several workers: Kuti and co-workers [19] have investigated the one-component scalar (j>^ model and 
the 0(4) symmetric scalar model and estimated a bound on the Higgs mass. Kroger and co-workers [20] have solved 
the Langevin equation on a momentum lattice for the scalar 4''^^i model and extracted critical behavior. Glueball 
scattering in compact QED2+1 (QC-D-like model) has been computed on a momentum lattice in Ref. [16]. Properties 
of nuclear matter have been computed by Brockmann and Frank [21]. Kogut and co-workers [22] have studied the 
phase diagram of quenched QED on a momentum lattice. Espriu [23] has studied the renormalization group flow by 
use a momentum lattice. Finally, Koutsoumbras [24] has computed the gluon propagator of flnite temperature QCD 
from a momentum lattice. Thus a momentum lattice regularization has proven useful when studying numerically 
physics near a critical point. 

When one treats gauge theories on a momentum lattice the following problem occurs: If one takes the gauge flelds 
Afj,{ki) as variables (so-called non-compact formulation), where ki denotes a momentum lattice, then the gauge action 
is not manifestly gauge invariant. As consequence one has observed non-local counter terms when computing from 
lattice perturbation theory the axial anomaly and the one-loop vacuum polarization. This has been seen by Karsten 
and Smit [37] by computing the triangle diagram using the SLAG derivative in the action and by Kroger and co- 
workers [38] using an action deflned on a momentum lattice with a momentum cut-off A. As Wilson has pointed out, 
it is desirable to conserve gauge symmetry manifestly in a regularized gauge theory. E.g., there is numerical evidence 
[39] that a lattice action being not manifestly gauge invariant yields no area law for the Wilson loop in pure SU (2) 
gauge theory. 

The space-time lattice Hamiltonian, corresponding to the Wilson action and being manifestly gauge invariant, has 
been constructed by Kogut and Susskind [40]. Here we are confronted with the following problem: How to introduce 
a momentum lattice as regulator while conserving manifestly gauge invariance? We suggest to do this as follows: We 
consider as variables closed Wilson-loops (for pure gauge theory without color charges). The Hilbert space is built 
from those loops. Gauge invariance corresponds to fulfllling Gauss' law 

G\ (!>>={). (51) 

For a flxed lattice site i, one has Gi = 'ij i i-^-i sum over generators of gauge transformations (where the 

temporal gauge is flxed). States corresponding to closed loops obey this law, while states corresponding open strings 
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do not obey it. Physical states are color singlet states, thus open string states are unphysical. Nevertheless we will 
make use of them as an intermediate step when constructing the Hilbert space of states obeying the Breit-condition. 

In order to introduce a regularization, we start from a conventional space-time lattice (regular, hypercube) with 
lattice spacing a. Then closed loops as well as open strings are defined as curves connecting adjacent lattice sites 
(straight line between neighboring lattice sites). E.g., a loop state is given by 

I (j>{xi) >=\ Ui_,{xi)U^{xi + afi) ■ ■ -UuiixN) >, (52) 
where x^^ + auj = xi. In order to introduce a momentum lattice we make a discrete Fourier transformation, 

I 4>i^i) >= a'^exp[-ixiki] \ <f){xi) >, (53) 

where each component of ki runs over the Brillouin zone —ir/a to +Tr/a. One can define the lattice momentum 
operator via the lattice translation Tfj_[a), which translates each configuration on the lattice by an increment a in 
the direction ji. It is given by 

T^(a) =exp[-iaP^]. (54) 
The eigenvalues of P are ki, which are the possible momenta of the loop state. 

In order to construct states with well defined momentum, obeying the Breit-condition, plus satisfying gauge in- 
variance, we suggest to proceed as follows: We construct a Hilbert space built from link states. By discrete Fourier 
transformation we associate a discrete lattice momentum to each link, say Ufi{ki). Then we construct multiple link 
states 

I ^{ki,---,kN) >=| U^{ki)UAk2)---U^{kN) > ■ (55) 
This state corresponds to momentum ktot = ki + ■ ■ ■ k^^ . Then we impose the Breit-condition (see Eq.(4)) 

ih - P/2f < {P/2)\ (56) 

where the parton momenta are given by the lattice momenta of the links. Thus, like in the scalar model, positivity 

\p\ 

of parton momenta in all three components and a given total momentum limits the total number of links to and 
thus give drastic bounds on the dimension of the effective Hilbert space. Eventually, we implement gauge symmetry 
by requiring Gauss' law, Eq.(51), to be respected. Thus the Breit-condition and Gauss' law define our basis of Hilbert 
states. 



VI. 5-MATRIX 



The Hamiltonian in the Breit-frame regularization has been shown above to be a suitable tool in the scalar model 
for computation of the mass spectrum and physics at the critical line as well as distribution functions. In this part 
we want to suggest that it is also a valuable tool for scattering phenomena and in particular for the non-perturbative 
computation of the S'-matrix. When considering non-perturbative computation of scattering observables, standard 
Euclidean lattice field theory is faced with the following problem: Scattering matrix elements are directly related to 
Minkowski n-point functions. On the lattice one can compute Euclidean n-point functions. In principle there is an 
analytic continuation between those two types of n-point functions. However, when the Euclidean n-point function is 
only known at some lattice points within the uncertainty of statistical errors, it is very difficult (almost impossible) 
to get reliable numerical results from an analytic continuation. A way out of this dilemma has been proposed by 
Liischer [29]. The idea is that continuum scattering phases can be extracted from the finite-size behavior of a mass 
spectrum on a finite lattice. This requires mass calculations via standard Euclidean lattice techniques, but requires 
quite precise data in order to resolve finite size effects. 

An alternative way to compute non-perturbatively an S'-matrix has been suggested by Kroger and co-workers [17]. 
The idea is the following. The S'-matrix, as has been introduced by Heisenberg [30] and M0ller [31], is defined as 
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5'=< V*"' I V-'^' >, 



(57) 



which is the probability amplitude to find an outgoing scattering state in an incoming scattering state. The scattering 
states are characterized by two conditions: (i) they are eigenstates of the Hamiltonian, and (ii) for t ±cxd they 
approach an asymptotic state. The asymptotic state describes two non-interacting particles (in the case of two-particle 
scattering). The so-called M0ller operator maps the asymptotic states | (f)""" > onto the scattering states | (f)'^^'> >. 

I ,^(±) >= fi(±) I >= s- lim exp[iHt] exp[-iH°t]<j,''' > . (58) 

i— )-=poo 

Those equations define scattering states and the S'-matrix. They can be taken over to quantum field theory with some 
care. 



A. Asymptotic one- and two-particle states 

One problem is the construction of an asymptotic one-particle state, asymptotic two-particle state, etc. In con- 
structive quantum field theory this is answered by Haag-Ruelle theory [32]. It tells how to construct asymptotic 
one-particle states by application of suitable local field operators on the physical vacuum. 

I 1 >phys= a^f) I >phys, (59) 

and a two-particle state is given by 

I 2 >phys= at(/i)at(/2) I >phy, . (60) 

Here is the creation operator of a one-particle state with wave function / created from the physical vaccum. The 
existence of such an operator has been proven by Haag and Ruelle [32]. An explicit form of this operator for the case 
of glueball states in pure gauge theory has been given by Liischer [33]. But Haag-Ruelle theory says nothing about 
how to find the physical vacuum. In the Hamiltonian approach in connection with the Breit-frame regularization, as 
advocated here, we avoid constructing the physical vacuum. Thus we proceed a route alternative to Haag-Ruelle's 
theory. We construct a one-particle state with momentum p directly by calculating an eigenvector of the regularized 
Hamiltonian H , 

H \p>= E{p} \p> . (61) 

The property of being a one-particle state is verified by computing its mass (see sect. 3). If, e.g., its mass is the lowest 
mass of the mass spectrum, the state | p > is a one-particle state. Let | / > denote such a one-particle state with 
momentum distribution given by a wavefunction /. In the language of Haag-Ruelle theory, the explicit construction 
of the state | / > means that we have found a creation operator (f) with 

|l>=At(/) |0>;,ee. (62) 

I.e., it creates a one-particle state from the vacuum of the regularized free Hamiltonian. There is a theorem by Haag 
[34] which says that in the continuum limit of relativistic quantum field theory, the physical Hilbert states of the 
interacting field theory (Hamiltonian) have nothing to do with those of the free field theory (free Hamiltonian). In 
particular there is no unitary transformation between the physical vacuum to the free vacuum. However, this theorem 
does not apply when we consider the regularized Held theory (Hamiltonian). Then there is a unitary transformation 
U , mapping the (regularized) free vacuum onto the (regularized) physical vacuum, 

I >;Zs= U I . (63) 

This relates the Haag-Ruelle creation operator (f) (of the regularized field theory) to the creation operator (f) 

by 

a\f) = UA\f)U-\ (64) 

Finally, using A^ (f) from Eq.(42) we can construct in equivalence to Eq.(40) asymptotic non-interacting two-particle 
states by 

|2>=At(/i)At(/2) |0>;,ee . (65) 
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B. M0ller wave operators and 5-matrix 



Let us denote by | (f)'" (pi, P2) > the asymptotic two-particle state, corresponding to two non-interacting particles 
with momentum pi and p2 respectively. Then the M0ller wave operator is given by 



i^s^LiT) >= n^^HT) I r'{pi,P2) >= exp[TiHT]exp[±i{E{p,) + E{p2))T] \ r'{pi 



P2) > ■ (66) 



Here E{p) denotes the energy-momentum dispersion relation of the one-particle state of mass m. H denotes the 
regularized Hamiltonian. The time parameter t, which goes to infinity in the continuum limit, has to be chosen to 
take a positive finite value T in the regularized theory (see below). In a similar way, one can construct the S'-matrix 

Sf.MT) =< ^f- I exp[i{E{p',) + E{p',))T] exp[-i2HT] exp[i{E{pi) + E{p2))T] \ > . (67) 

From the numerical point of view, the computation of the S'-matrix element proceeds most simply by diagonalizing 
the regularized Hamiltonian 

H \ T]^ >= \ T]^ >, v = 1,2, - ■ ■ 

T]^ > exp[iej.T] < r]^ \ ■ (68) 

How should one choose the scattering time parameter T? When applying this time-dependent Hamiltonian method 
to non-relativistic quantum mechanics as well as to field theoretic models [17], the following general observations have 
emerged from numerical calculations: The matrix element Sfi^in{t) considered as a function oft has the following 
behavior. At t = it takes the value < (fij' \ (f)"^ > (in the case of elastic scattering). When increasing t it deviates 
from the starting value and eventually reaches a plateau region. When further increasing t, it leaves the plateau region 
and after a while shows an (irregular) oscillatory behavior. The plateau region is the region of physical interest. Its 
existence can be shown analytically for non-relativistic potential scattering (see Ref. [17] and references therein). 
The location and size of this plateau region depends on the model and dimension. In particular, it depends on the 
dimension of the regularized Hamiltonian. When increasing this dimension, i.e., when exploring a larger Hilbert space, 
the size of the plateau region becomes larger. In the continuum limit, when the S'-matrix converges, the size of the 
plateau should become infinitely large. The time parameter T should be chosen from this plateau region, either by 
determining where the matrix element Sji^in (t) has the least variation in t, or by the following criterion of conservation 
of energy: In the continuum limit, energy conservation in a scattering reaction means that 

<^(±) >=E„,, (69) 

where Eas denotes the energy of the asymptotic non-interacting two-particle state. Thus we define the function 



AE{t) 



< n^^Ht)ct,-^p,,p2) I H I fi(±)(t),/)«^(pi,p2) > -Eas /Eas, (70) 



where Sl'-^-' is given by Eq.(58). This function is a measure of violation of energy conservation in a scattering reaction 
computed with the regularized Hamiltonian at some finite time t. In the continuum limit this should be zero. Thus 
we can choose the time paramer T such that AE{t) has a minimum. Numerical experience has shown that the value 
of T once determined as position of minimal variation of the S'-matrix element and secondly determined as position 
of the minimum of AE are quite close together, which is an indication of consistence. 



In order to get the physical S'-matrix one has to carry out renormalization and take into account the vacuum 
structure. Renormalization means that firstly one has to determine the counter terms in the Hamiltonian. E.g., 
for the scalar (j>^ model, one has to renormalize the wavefunction, the mass and the coupling constant. Then one 
computes physical observables like, e.g., masses or scattering cross sections and tunes the bare parameters of the 
model, such that the physical observables remain fixed. Finally, the vacuum structure needs some careful treatment. 
The computation of the S'-matrix, as described above yields the full S'-matrix, which includes the connected part 
(which is the piece observed in scattering experiments) but also all disconnected parts. The factorization of n-point 
Greens functions into connected pieces, is knows as vacuum structure [35]. This allows to extract the connected part 
of the S'-matrix. 
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The time-dependent Hamiltonian method as described above, but instead of using the Breit-frame regularization by 
use of the rest-frame regularization, has been applied to glueball scattering in compact U{1) gauge theory (compact 
QED) in 2 -|- 1 dimensions [16]. At the end of this section we want to adress the question: What advantage does 
it bring to use the Breit-frame regularization for scattering calculations in the time-dependent Hamiltonian formu- 
lation? Firstly, as mentioned above, the Breit-frame regularization avoids the calculation of the vacuum state when 
constructing asymptotic non-interacting two-particle states. Secondly, this regularization reduces the number of ef- 
fective degrees of freedom by the same mechanism as was shown to be useful for the calculation of the mass spectrum. 
However, one must pay attention to the following limitation: Because we take into account parton momenta inside 
the sphere given by the Breit condition Eq.(4), the momenta of the asymptotic particles, i.e., p'l, p'2, Pi, P2 should lie 
well inside the Breit sphere. This limits the scattering reactions which can be treated. E.g., head-on collisions are 
not included. However, this constraint is not very stringent, because a suitable Lorentz-boost can be applied to map 
the momenta into the Breit sphere. 



VII. FINITE DENSITY THERMODYNAMICS 



The computation of thermodynamic observables at finite temperature and finite density is an important problem 
for the physics of neutron stars, high energy heavy ion collisions, and for the question of phase transition from the 
hadronic phase to quark-gluon plasma in QCD. However, when treating finite temperature QCD in the standard 
Lagrangian lattice approach, there is an well known problem when a non-zero chemical potential is included to describe 
finite density effects. Then the fermionic determinant becomes complex yielding a complex lattice action. This has 
led to great difficulties when solving the model numerically via Monte Carlo methods [36]. In order to study the 
infrared dynamics of Yang-Mills and Yang-Mills-Higgs theories at finite temperature, which can not be addressed by 
Euclidean methods, Moore [41] has suggested an improved Hamiltonian for Minkowski Yang-Mills theory. 

In this section we want to discuss how finite temperature and finite density thermodynamics can be treated in a 
Hamiltonian formulation with the Breit frame regularization. The point is that the Hamiltonian formulation allows 
to treat also non-Hermitian Hamiltonians (complex actions). Consider the following partition function 

Z = Trexp[--^iHo+iiN)], (71) 

where Hq stands for a Hermitian Hamiltonian, /j, denotes the chemical potential and N stands for a particle number 
operator. Let us suppose now, for the sake of the argument that the term /^N would be non-Hermitian. What is 
then the advantage to use a Hamiltonian formulation? In a Hamiltonian formulation this partition function can be 
computed non-perturbatively via diagonalization of Hq + jJ-N in the same way as exp[iHt\ has been computed when 
calculating the S'-matrix (sect. 6). 

What is the advantage to use the Breit-frame regularization? Let us consider the following scenario: One wishes 
to study hot nuclear matter, respectively quark-gluon plasma in a state with total momentum P z/z 0. In order to 
investigate this experimentally, one can perform a high energy heavy ion collision experiments with large momentum 
P. However, it is stil an open question if thermodynamical equilibrium can be reached during the short time of 
collision and before decaying of fragments, which would justify the use of the Boltzmann -Gibbs partition function. 
Putting aside for the moment the question of experimental realization, it is nevertheless physically interesting to ask 
the following question: What are the properties of matter at finite temperature and finite density at thermodynamical 
equilibrium in a sector of momentum P z/z 0. Thus we consider the partition function at momentum P, 



ZIP) = Tr 



krI 



(72) 

P 



The evaluation of this function now can be done in the Breit-frame regularization, which by the same reason as in 
the computation of structure functions would reduce the effective number of degrees of freedom, i.e., the dimension 
of the effective Hilbert space. 



VIII. SUMMARY 



In conclusion, we have suggested a Hamiltonian method and a momentum regularization corresponding to the 
Breit-frame. We have shown that this method allows to extract continuum physics by presenting numerical results 
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for the 4't+i theory in the symmetric phase close to the critical line. We find close agreement with the solution of the 
renormalization group equations by Liischer and Weisz. We have seen scaling behaviour of several low-lying masses 
near the critical point. Using the Breit-frame, we have computed analytically for DIS in QCD the relation between 
the hadronic tensor, the structure functions and the quark distribution functions. In the Bjorken limes we find the 
conventional relations between Fi, F2, gi and the quark distribution functions. We have presented numerical results 
for parton distribution functions for the <f)'^ and <f)^ model. The example of the <f)^ model demonstrates how a peak at 
small xb can be produced, with a behavior similar to Altarelli-Parisi behavior in QCD. We have proposed how the 
Breit-frame regularization can be applied to gauge theories, while keeping gauge symmetry manifestly conserved. We 
have suggested that this regularization might be useful also for the computation of scattering reactions (S'-matrix), 
as well as finite temperature and finite density thermodynamics. We are optimistic that the method can be applied 
to compute numerically structure functions in QCD. Work is in progress. 
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Fig.l The ground state mass niR in lattice units (a = 1) versus k for A = 0.00345739 (A = 0.01 in Ref. [2]), A and 
K are given by Eq.(6). The points correspond to results of Ref. [2]. Our results correspond to A/ Ap = 3 (dashed line) 
and A/ Ap = 4 (solid line). 

Fig. 2 The lowest lying mass spectrum versus k. The ground state mass is set to one, A as in Fig.[l]. A and k are 
given by Eq.(6). A/ Ap = 4. 

Fig. 3 The distribution function /{xb) of (j>^_^_-^ versus the momentum fraction xb- A = 0.00345739 (as in Fig.[l]), 
A/Ap = 4. 

Fig. 4 The distribution function f{xB,go) of <f)i_^_i versus the momentum fraction xb and the coupling constant go. 
The bare mass mo has been set to mo = 3Ak. A/ Ap =11. 
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IX. TABLE CAPTION 



A 


0.0005 


0.001 


0.005 


0.01 


0.05 


0.1 


'^crit 


0.125101 


0.125202 


0.125991 


0.126968 


0.132368 


0.13601 


a 


0.99997 


0.99993 


0.99972 


0.9993 


1.0073 


1.0275 



The critical points Kcrit versus A. A and k are given by Eq.(6). k^^^. is taken from Ref. [2]. a := K^rft/'^crit denotes 
the ratio between the results of this work and Ref. [2]. In this work, Kcrit has been determined under the condition 
that the renormalized mass niR becomes imaginary. A/ Ap = 4. 
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